Current-carrying molecules: a real space picture 
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An approach is presented to calculate characteristic current vs voltage curves for isolated molecules 
without explicit description of leads. The Hamiltonian for current-carrying molecules is defined by 
making resort to Lagrange multipliers, while the potential drop needed to sustain the current is 
calculated from the dissipated electrical work. Continuity constraints for steady-state DC current 
result in non-linear potential profiles across the molecule leading, in the adopted real-space picture, 
to a suggestive analogy between the molecule and an electrical circuit. 
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Experiments on single-molecule junctions are challeng- 
ing, mainly due to the need of contacting a microscopic 
object, the molecule, with macroscopic leads |il(, .2', 3, 4]. 
Theoretical modeling of molecular junctions is difficult 
d, 0, 0, 13 , and again the description of contacts repre- 
sents a delicate problem. To attack the complex prob- 
lem of conduction through a molecular junction a strat- 
egy is emerging |^ ^| that focuses attention on isolated 
molecules and describes the intrinsic molecular conduc- 
tivity in the absence of electrodes. At variance with 
common approaches that impose a potential bias to 
the electrodes and then calculate the resulting current 

0, 0, , a steady-state DC current is forced in the iso- 
lated molecule by making resort to a Lagrange-multiplier 
technique 0, or by drawing a magnetic flux through 
the molecule |0. Whereas the strategy is promising, 
two main problems remain to be solved: (1) the calcu- 
lation of the potential drop needed to sustain the cur- 
rent, and (2) the definition of the potential profile in the 
molecule. Here I demonstrate that the Joule law can be 
used to calculate the potential drop from the electrical 
power dissipated on the molecule. Moreover, continuity 
constraints for steady-state DC current are implemented 
in polyatomic molecules in terms of multiple Lagrange 
multipliers that yield to non-linear potential profiles in 
the molecule. Finally, in the adopted real-space picture, 
the current flows through chemical bonds rather than 
through energy levels, leading to a suggestive description 
of the molecule as an electrical circuit with resistances 
associated to chemical bonds. 

To start with consider a diatomic Hubbard molecule, 
whose Hamiltonian Hq is defined by U, t, and the differ- 
ence of on-site energies: 2 A — €2 — £i- Following Kosov 
, a current is forced through the molecule by introduc- 
ing a Lagrange multiplier. A, as follows: 



(1) 



where cj ^ creates an electron with spin a on the i-site, 

and j — —it'^^{c\^C2a — H .c.) measures the current flow- 
ing through the bond. Here and in the following h and 
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the electronic charge are set to 1, and t is taken as the 
energy unit. The ground state of H{X), \G{X)), carries a 
finite current, J — {G{X)\j\G{X)), and the Lagrange mul- 
tiplier, A, is fixed by imposing a predefined J 0. Other 
molecular properties can be calculated as well, and their 
dependence on J can be investigated 9]. Just as an ex- 
ample, the bond-order decreases with J, and, in systems 
with inequivalent sites, the on-site charge distribution is 
equalized by the current flow. These are interesting infor- 
mations, but characteristic J{V) curves are still needed. 

The Lagrange multiplier. A, has the dimensions and the 
meaning of a magnetic flux drawn across the molecule to 
generate a spatially uniform electric field [lol IT]| , E cx 
LuX, where iv is the field frequency In the limit of 

static fields, uj 0, both E and V vanish, suggesting that 
a finite current flows in the molecule at zero bias. This 
contrasts sharply with the fundamental relation between 
charge transport and energy dissipation ^ ^j, |0: a 
finite V is needed to sustain a current due to dissipative 
phenomena occurring in the conductor. Specifically, the 
Joule law relates the potential drop in a conductor to 
the electrical power spent on the system to sustain the 
current, W — VJ. Since J is known, V can be obtained 
from a calculation of the dissipated power. 

Dissipation is conveniently described in the density 
matrix formalism using as a basis of the eigenstates \k) of 
H{0) n^H^ . The equilibrium density matrix, ctq, is a di- 
agonal matrix whose elements are fixed by the Boltzmann 
distribution. On the same basis (t(A) is a non-diagonal 
matrix corresponding to a non-equilibrium state whose 
dynamics is governed by: a = —j:[H, <t] + &r, where an 
accounts for relaxation phenomena, as due to all degrees 
of freedom not explicitly described by H (e.g. molecu- 
lar vibrations, or environmental degrees of freedom also 
including leads) Hill 111. Diagonal elements of ctr de- 
scribe depopulation and are associated with energy dissi- 
pation. As for depopulation I adopt a simple phenomeno- 
logical model with {&R)kk = -J^mlkmO-mm, where jkm 
measures the probability of the transition from k to m 
[1^ IT3 |. For the sake of simplicity I will consider the 
low-temperature limit, with jkm = 7 for k > m and 
7fcm = otherwise. The dynamics of off-diagonal ele- 
ments is governed by depopulation and dephasing effects: 

(O'fi)fcm = -rfcmCTfem, with Tkm = {jkk + "/mm)/'^ + 7fcm ' 
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FIG. 1: Characteristic J{V) curves for a two-site system with 
different A and U. The depopulation rate 7 is set to 0.2, 
dephasing is neglected. 



where jkk = J2'mlrnk, and 7^^ describes dephasing i.e. 
the loss of coherence due to elastic scattering |0, QJ] . 

The energy dissipated by the system, Tr{&rH), has 
two contributions: the first one, Wd = X]fc(^Jj)fefci 
ways negative and measures the energy that the system 
dissipates to the bath as the current flows. This term 
is governed by depopulation, whereas dephasing plays 
no role. The second term, W — —XTr{(7iij), measures 
the electric work done on the system to sustain the cur- 
rent: it is this term that enters the Joule law. In non- 
degenerate systems j is an off-diagonal operator, so that 
only off-diagonal elements of &r enter the expression for 
W. Both depopulation and dephasing then contribute 
to y, and hence to the molecular resistance. This is in 
line with the observation that a current flowing through 
a molecule implies an organized motion of electrons along 
a specific direction f 2,'l5]. Therefore any mechanism 
of scattering, either anelastic, as described by depopu- 
lation, or elastic, as described by dephasing, contributes 
to the electrical resistance 0|. The unbalance be- 
tween Wd and W is always positive: the molecule heats 
as current flows. Efficient heat dissipation is fundamen- 
tal to reach a steady-state regime and to avoid molecular 
decomposition ,,12.] . 

Fig. 1 shows the characteristic curves calculated for 
a diatomic molecule with 7 = 0.2 and 7^^ = 0. In the 
left panel results are shown for the symmetric, A = 0, 
system. As expected, electronic correlations decrease the 
conductivity. The results in the right panel for an asym- 
metric system (A ^ 0) show instead an increase of the 
low-voltage conductivity with increasing U. This inter- 
esting result is related to the minimum excitation gap, 
and hence the maximum conductance, of the system with 
U = 2 A. The asymmetric diatomic molecule represents a 
minimal model for the Aviram and Ratner rectifier |l6j |. 
however the characteristic curves in the right panel of 
Fig. ^are symmetric, and do not support rectification. 
In agreement with recent results, rectification in asym- 
metric molecules is most probably due to contacts ■^j, 
or to the coupling between electrons and vibrational or 
conformational degrees degrees of freedom p^ . 

Before attacking the more complex problem of poly- 
atomic molecules, it is important to compare the results 



obtained so far with well known results for the optical 
conductivity ^ ^km = F, as it occurs, e.g., for sys- 
tems with large inhonrogeneous broadening (the coherent 
conductance limit 0,0]), the expression for the poten- 
tial drop is very simple: V = XT. Then, a perturbative 
expansion of J leads to the following expression for the 
zero-bias conductivity, Qq: 
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(2) 



where g is the gs of H{0), with energy Eg, and the sum 
runs on all excited states. This expression for the DC 
conductivity coincides with the zero-frequency limit of 
the optical conductivity jllj , provided that the frequency, 
uj, appearing in the denominator of the expression for 
the optical conductivity in Ref. is substituted by 

Lu — iV. Introducing a complex frequency to account for 
relaxation is a standard procedure in spectroscopy |l4j . 
leading to similar effects as the introduction of an ex- 
ponential switching on of the electromagnetic field '■ 
both phenomena account for the loss of coherence of elec- 
trons driven by an EM field and properly suppress the 
divergence of the optical conductivity due to the build-up 
of the phase of electrons driven by a static field. 

The connection between DC and optical conductiv- 
ity breaks down in polyatomic molecules. To keep 
the discussion simple, I will focus attention on linear 
Hubbard chains. The optical conductivity of Hubbard 
chains was discussed based on the current operator J — 
'~ieY^^ti{c\ ^Ci+i,cr — H.c.)/h However, this opera- 
tor measures the average total current and does not ap- 
ply to DC currents. Specifically, if a term — A J is added 
to the molecular Hamiltonian, a finite average current 
is forced through the molecule 0, but this current does 
not satisfy basic continuity constraints for steady-state 
DC current. In fact, to sustain a steady-state DC cur- 
rent one must avoid the build up of electrical charge at 
atomic sites. Specifically, in linear molecules the conti- 
nuity constraint imposes that exactly the same amount 
of current flows through each bond in the molecule. To 
impose this constraint the current on each single bond 
must be under control and a Lagrange multiplier must 
be introduced for each bond, as follows: 



Xijz 



(3) 



where ji — ~iti{(^^ ^Cij^i^cr — H.c.)/h, and the A^'s are 
fixed by imposing ji — {G\ji\G) = J independent on i. 
As before, the electrical work done on the molecule is: 



(4) 



that naturally separates into contributions, Wi, rele- 
vant to each bond. The total potential drop across the 
molecule, V = W/J, is then the sum of the potential 
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FIG. 2: Left panel; characteristic curve for the three-site 
Hubbard chain sketched in the figure, with 3 electrons, U = 4, 
constant on-site energies and ti — 1.2, t2 — 0.8. Dephas- 
ing is set to zero, and depopulation rate is 7 = 0.1. Right 
panel: total resistance (R) and bond-resistances (_Ri and -R2). 
For bond-resistances continuous and dashed lines show results 
obtained by allowing the current to flow through the whole 
molecule, and through a single bond, respectively. 



drops across each bond, Vi = Wi/J, leading in gen- 
eral to non-linear potential profiles. Of course, in the 
adopted real-space picture the potential profile can only 
be calculated at atomic positions, and no information 
can be obtained on the potential profile inside each bond. 
Therefore, instead of showing the potential profile along 
the molecule, I prefer to convey the same information in 
terms of bond-resistances, defined as: Ri = [dJ/dVi)~^. 
Fig. 13 shows the behavior of a 3-site chain with three 
electrons, U = 4, equal on-site energies and different t: 
ti = 1.2, and t2 — 0.8. The left panel shows the char- 
acteristic J{V) curve, and continuous lines in the right 
panel report the total resistance R, and the two bond 
resistances, Ri and i?2. The molecular resistance varies 
with the applied voltage and, as expected, the resistance 
of the weaker bond is higher than the resistance of the 
stronger bond. 

Dimensionless resistances in the figure are in units with 
h/e^ = {2TrgQ)~^ = 1, where go is the quantum of conduc- 
tance, that, in standard approaches to molecular junc- 
tions jy| represents the maximum conductance associated 
with a discrete molecular level. This well known result 
is related to the inhoniogeneous broadening of molecular 
energy levels as due to their interaction with the elec- 
trodes |6[. Of course there is no intrinsic limit to the 
conductivity in the model for isolated molecules discussed 
here. 

As a direct consequence of the continuity constraint, 
the total resistance R — {dJ/dV)~^ is the sum of the 
two bond-resistances, leading to a suggestive description 
of the molecule as an electrical circuit, with resistances 
associated with chemical bonds joint in series at atomic 
sites. Whereas this picture is useful, the concept of bond- 
resistance should be considered with care in molecular 
circuits. At variance with standard conductors, in fact, 
the resistance of the bonds depends not only on the cir- 
cuit (the molecule) they are inserted in, but also on the 
way the resistance is measured. Dashed lines in Fig. |5| 



FIG. 3: The same as in figH but with 7 = and 7fc,„ = 0.2. 

show the bond-resistances calculated by forcing the cur- 
rent through specific bonds (i.e. by setting a single / 
in Eq. and these differ from the bond-resistances 
calculated when the whole molecule carries the current 
(continuous lines). 

The situation becomes somewhat simpler in the coher- 
ent conductance limit, Tkm = F. Perturbative arguments 
can be used to demonstrate that at zero bias the bond 
resistances calculated for the current flowing through the 
whole molecule or through a single bond do coincide. 
This additive result for the molecular resistance in the 
coherent transport limit is in line with the observation of 
transmission rates inversely prop ortional to the molecu- 
lar length in the same limit However, as shown in 
Fig. 121 this simple Ohmic behavior breaks down quickly 
at finite bias. 

The introduction of as many Lagrange multipliers 
as many current-channels (bonds) are present in the 
molecule accounts for a non-linear potential profile 
through the molecule, i.e. for a non-uniform electric field. 
Accounting for a single Lagrange multiplier coupled to 
the total current operator is equivalent to draw a mag- 
netic flux through the molecule as to generate a spatially 
homogeneous electric field |lOl ITll Il8j| . a poor approx- 
imation for DC conductivity in extended (polyatomic) 
molecules. Just as an example, for a 4-site chain with 
the same ti — 1 on each bond, the zero-bias resistance of 
the central bond exceeds that of the lateral bonds with 
R2/R1 ranging from 10 to 2 as J7 increases from to 
4. Bonds with the same t have different resistances due 
to their different bond-orders, and, in agreement with 
recent results 0, 1131 , this demonstrates nicely the need 
of accounting for non-uniform electric fields in extended 
molecules, even for very idealized molecular structures. 

It is of course possible to discuss more complex molec- 
ular models. As an interesting example, a nearest- 
neighbor hopping t' is added to the Hamiltonian for the 
three site molecule discussed above. This opens a new 
channel for electrical transport, and a term —X'j' adds 
to the Hamiltonian with j' — —it' ^^{c\^C3a — H.c ). 
As before, continuity imposes ji — j2, as to avoid build- 
ing up of charge at the central site. The total current 
is J = ji -f j' and, of course, no continuity constraint 
is given on j' . However the potential drop across the 
molecule, i.e. the potential drop measured at sites 1 and 
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FIG. 4: Left panel: characteristic curves of the three site 
chain sketched in the figure, with three electrons, constant 
on-site energies, U — t\ — t2 — \- and '^um = and — 
0.2. Continuous and dotted lines refer to t' = 0.4 and 0, 
respectively. Right panel: molecular and bond-resistances for 
the chain with t' = 0.4. The dotted line shows the total 
resistance for the chain with t' = 0. 



3 must be uniquely defined. Therefore one must tune 
Ai, A2 and A' as to satisfy j\ — ji^ while satisfying the 
condition: V\ ^- Vi = V — V, with Vi — Wi/ji and 
V' = W /j'. Imposing a constraint on the potentials is a 
tricky affair, that becomes trivial when Tkm — T- In that 
case in fact Vi = TXi and V = FA', so that the constraint 
on the potentials immediately translates into a constraint 
on Lagrange multipliers. Fig. 0] shows some results ob- 
tained in this limit for a system with ti = t2 = 1, t' = 0.4. 
In spite of the fairly large t' value, the contribution to 
the current from the bridge-channel is small, mainly due 
to the small bond-order for next-nearest neighbor sites. 
Once again the physical constraints imposed to the cur- 
rents and to the potentials lead to standard combination 
rules for bond-resistances with 1/R= 1/ R' +l/{Ri+R2). 
As for the DC conductivity is concerned, the molecule be- 
haves as an electrical circuit with two resistances, i?i and 
i?2 in series bridged by a parallel resistance, i?'. 

Applying the proposed approach to complex molec- 



ular structures and/or to molecules described by accu- 
rate quantum chemical Hamiltonians is non-trivial due 
to the appearance in the Hamiltonian of as many La- 
grange multipliers as many current channels are consid- 
ered, and due to the large number of constraints to be 
implemented. Instead, at least for small molecules, the 
approach can be fairly easily extended to account for vi- 
brational degrees of freedom. Non-adiabatic calculations 
are currently in progress to describe electrical conduction 
through a diatomic molecule in the presence of Holstein 
and Peierls electron-phonon coupling. More refined mod- 
els for the relaxation dynamics can also be implemented 
|19| , whereas the introduction of spin-orbit coupling can 
lead to a model for spintronics. 

In conclusion, this paper presents an approach to the 
calculation of characteristic current/ voltage curves for 
isolated molecules in the absence of contacts. While 
hindering the direct comparison with experimental data, 
this allows the definition of the molecular conductivity as 
an intrinsic molecular property. Even more important, 
a paradigm is defined for imposing a steady-state DC 
current through the molecule, while extracting the volt- 
age drop across the molecule from the energy dissipation. 
The careful implementation of continuity constraints for 
steady-state DC current leads to the definition of the po- 
tential profile through the molecule, that, in the adopted 
real-space description, quite naturally results in the con- 
cept of bond-resistances, in a suggestive description of 
the molecule as an electrical circuit with current flowing 
through chemical bonds. 
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